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PREFACE 


The  angular  correlation  of  annihilation  radiation  by 
slit  geometry  can  be  used  to  obtain  the  momentum  distribu¬ 
tion  of  annihilating  electron-positron  pairs  in  matter. 

The  purpose  of  this  thesis  is  to  describe,  evaluate  and 
compare  two  types  of  apparatus  employing  rectangular  slits 
and  two  methods  of  data  analysis  that  can  be  used  in  this 
work. 
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CHAPTER  I 


DESCRIPTION  OF  THE  TYPES  OF  APPARATUS 

When  a  positron-electron  pair  annihilates  by  two 
photon  annihilation,  the  photons  will  be  emitted  in  ex¬ 
actly  opposite  directions  if  the  pair  is  initially  at 
rest.  A  deviation  of  the  angle  between  the  emitted  pho¬ 
tons  from  180°  indicates  an  initial  component  of  momen- 
turn.  This  deviation  from  180°  has  rarely  been  found  to 
exceed  15  milliradians.  Angular  correlation  of  these 
photons  over  this  range  yields  information  about  the 
momentum  distribution  of  the  annihilating  pairs. 

The  experimental  configuration  that  employs  rec¬ 
tangular  detectors  is  referred  to  as  slit  geometry. 

Two  types  of  apparatus  employing  slit  geometry  have 
been  used.  One  type^»^  consists  of  two  very  long  and 
narrow  slits  with  a  single  coincidence  channel  while 
the  other  type^’^  consists  of  eight  narrow  slits 


Berko,  S.,  and  Plaskett,  J.  S, ,  Phys.  Rev.  112. 
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(four  at  each  end)  with  sixteen  coincidence  channels. 

It  was  originally  planned  to  construct  a  multiple  coin¬ 
cidence  channel  apparatus  here,  but  this  plan  was  later 
changed  because  of  the  reasons  that  will  be  discussed 
in  this  thesis.  Some  work  was  done  on  this  type  of 
apparatus,  however,  and  Figure  1  shows  a  sketch  of  the 
original  plan  modeled  after  Trumpy's^  apparatus.  This 
was  later  modified  for  single  coincidence  channel  work, 
and  a  sketch  of  this  apparatus  is  shown  in  Figure  2. 

The  major  difference  between  the  two  instruments 
is  the  slit  and  photomultiplier  tube  arrangement.  In¬ 
strument  1  (Figure  1)  was  to  have  used  eight  RCA  6342A 
photomultiplier  tubes  with  a  1%  inch  diameter  by  \  inch 
thick  Nal(Tl)  crystal  mounted  directly  on  the  photoca¬ 
thode  of  each.  Four  of  these  were  to  be  mounted  at  each 
end  of  the  instrument  with  lead  shielding  placed  in  such 
a  way  that  only  a  lower  slice  of  each  crystal,  the  thick¬ 
ness  of  which  would  subtend  an  angle  of  0.4  milliradian 
at  the  target,  could  detect  radiation.  The  crystal 
slices  at  the  stationary  end  were  to  have  been  0.6  mil¬ 
liradian  apart,  and  those  at  the  movable  end  were  to 
have  been  2.4  milliradians  apart  making  it  possible  to 
simultaneously  measure  16  values  of  0  covering  9.6  mil¬ 
liradians  in  steps  of  0.6  milliradian.  Instrument  2 
(Figure  2)  uses  two  RCA  6342A  photomultiplier  tubes, 
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Figure  1 

Side  View  of  16  Coincidence 
Channel  Apparatus 
(Not  to  Scale) 


SIDE  VIEW  OF  16  COINCIDENCE  CHANNEL  APPARATUS 


Figure  2 

Top  View  of  Single  Coincidence 
Channel  Apparatus 
(Not  to  Scale) 


TOP  VIEW  OF  SINGLE  COINCIDENCE  CHANNEL  APPARATUS 
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each  looking  into  the  end  of  a  1  inch  diameter  by  8  inch 
long  Nal(Tl)  crystal.  One  tube  and  crystal  are  placed 
in  lead  shielding  at  each  end  of  the  instrument  in  such 
a  manner  that  an  8  inch  by  .032  inch  portion  of  the  crys¬ 
tals  is  visible  to  the  targets.  The  result  is  an  instru¬ 
ment  that  employs  one  very  long  and  very  narrow  slit  at 
each  end. 

Three  32  inch  high  cinder  block  piers  were  constructed 
along  a  straight  line  2  meters  apart.  The  center  pier 
supports  the  source  and  target  holder,  and  the  other  two 
piers  support  the  detectors.  An  86  inch  long  platform 
made  of  aluminum  was  constructed  so  that  one  end  pivots 
on  a  support  at  the  central  pier  and  the  other  end  tra¬ 
vels  up  and  down  on  a  precision  screw  mounted  on  one  of 
the  outer  piers.  The  pitch  of  the  screw  is  such  that 
one  full  turn  will  increase  or  decrease  the  angle  that 
the  platform  makes  with  the  horizontal  by  0.6  millira- 
dian.  The  screw  is  linked  to  a  four-pole  Geneva  Movement 
by  one  to  four  gear  ratio  so  that  one  operation  of  the 
Geneva  Movement  will  turn  the  screw  through  360°.  The 
Geneva  Movement  is  powered  by  a  reversible  electric  mo¬ 
tor  that  can  be  programmed  to  turn  the  screw  whenever  a 
change  in  the  angle  between  the  slits  is  desired.  A 
sketch  of  this  mechanism  is  shown  in  Figure  3. 

The  source  and  target  holder  (see  Figure  4)  con¬ 
sists  of  a  lead  cylinder  10  inches  in  diameter  by  10 
inches  long.  Suitable  means  for  adjusting  the  orientation  of 


Figure  3 

Sketch  of  Screw  and  Geneva 
Movement  Mechanism 
(Not  to  Scale) 


SKETCH  OF  SCREW  AND  GENEVA 
MOVEMENT  MECHANISM 

(NOT  TO  SCALE) 


Figure  4 

Sketch  of  Source  and  Target  Holder 
(Not  to  Scale) 
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the  holder  are  provided  so  that  its  axis  can  be  aligned 
with  the  axis  of  the  instrument.  The  source  is  mounted 
on  the  end  of  a  brass  cylinder  (S)  that  is  inserted  from 
the  top.  The  depth  of  penetration  is  controlled  by  two 
brass  screws  (Sc^).  The  target  is  mounted  on  a  circular 
brass  plate  (Ta)  that  is  connected  to  the  face  of  a 
brass  cylinder  (T)  by  three  equally  spaced  spring-loaded 
screws.  This  arrangement  allows  an  orientation  of  the 
target  that  is  independent  of  the  orientation  of  the 
source  and  target  holder.  The  target  cylinder  fits  into 
the  source  and  target  holder  from  the  bottom,  and  the 
brass  cylinder  rests  upon  a  screw  (Sc^)  that  controls 
the  depth  of  penetration.  The  collimators  are  two  lead- 
filled  brass  cylinders  (C)  having  rectangular  slits  (SI) 
machined  in  the  lead.  These  cylinders  fit  into  each 
end  of  the  source  and  target  holder.  O-ring  grooves 
are  provided  wherever  necessary  so  that  the  interior  of 
the  source  and  the  target  holder  can  be  evacuated  if 
later  desired. 

Solid  state  electronics  is  being  used  wherever 
possible.  The  output  of  each  photomultiplier  tube  goes 
to  a  preamplifier  and  pulse  shaper.  The  outputs  from  the 
pulse  shapers  go  to  the  coincidence  circuit,  and  the 
coincidence  counts  are  scaled  down  by  a  factor  of  four. 
The  output  of  the  scale  of  four  is  then  directed  to  the 
registers  through  a  driver  and  switch.  Schematics^ 
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The  development  of  the  solid  state  electronics 
used  in  this  apparatus  was  done  by  Mr.  Louis  Deiterman. 


17 


of  the  preamplifier,  pulse  shaper,  and  coincidence  cir¬ 
cuit  are  indicated  in  Figure  5,  6  and  7  respectively. 

The  original  plan  (16  coincidence  channels)  called  for 
a  complex  switch  to  direct  the  output  of  each  coincidence 
circuit  to  the  correct  sequence  of  registers.  This  switch 
was  designed  and  is  shown  in  Figures  8  and  9  in  case  it 
might  be  of  future  use. 


Figure  5 

Schematic  of  Preamplifier 
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a. 


SCHEMATIC  OF  PREAMPLIFIER 


Figure  6 

Schematic  of  Pulse  Shaper 


SCHEMATIC  OF  PULSE  SHAPER 


Figure  7 

Schematic  of  Coincidence  Circuit 


Figure  8 

Schematic  of  Switch 
(T-Bar  Relays  and  Code) 


Figure  9 

Schematic  of  Switch 

(Control  System  for  Coding  T-Bar  Relays) 


Coincidence 


SCHEMATIC  OF  SWITCH 
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CHAPTER  II 


RESOLUTION  AND  DATA  ANALYSIS 

The  momentum  distribution  of  annihilation  photons 
can  be  obtained  from  the  data  in  two  ways.  The  first 
method  to  be  discussed  (Method  A)  is  to  determine  the 
resolution  of  the  instrument  in  the  z  direction  and  from 
the  data  corrected  in  this  manner  to  calculate  the  dis¬ 
tribution  in  the  z  component  of  the  momentum.  The  second 
method  to  be  discussed  (Method  B)  is  to  determine  the  Q 
resolution  of  the  instrument  (where  -  Q  is  the  angle 
between  the  photons)  and  from  the  data  corrected  in  this 
manner  to  calculate  the  distribution  in  the  p  -component 

(where  p  =  -/  p2  + p*  )  of  the  momentum.  The  results 
p  x  i 

obtained  from  each  method  allow  calculation  of  the  total 
momentum  distribution  of  annihilation  photons. 

The  following  discussion  will  refer  to  a  coordinate 
system  that  is  centered  at  the  target  with  the  z-axis 
vertical  and  the  y-axis  along  the  horizontal  line  defined 
by  the  source  and  the  stationary  slit. 
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A.  Method  A 

This  method  has  been  used  by  several  authors^ >8  and 
will  be  discussed  in  some  detail.  The  finite  resolution 
of  the  apparatus  is  a  result  of  the  finite  dimensions 
of  the  target  and  of  the  slits.  Since  the  z  component 
of  the  momentum  is  being  measured,  it  is  desirable  to 
have  as  narrow  a  resolution  function  in  the  z  direction 
as  possible.  The  calculation  of  the  resolution  function 
will  show  that  this  requirement  makes  it  necessary  to 
use  a  target  and  slits  that  are  very  narrow  in  the  z 
direction.  The  stationary  slit  will  be  referred  to  as 
Slit  A  and  the  movable  slit  will  be  referred  to  as  Slit  B. 
Each  slit  has  the  dimensions  of  2bQ  in  the  z  direction 
and  2aQ  in  the  x  direction  and  is  a  distance  rQ  from  the 
target. 

If  an  annihilation  pair  at  the  origin  have  a  mo¬ 
mentum  component  in  the  z  direction  of  p2  and  if  one 
photon  passes  along  the  y-axis  through  the  center  of 
Slit  A  (see  Figure  10),  the  other  photon  must  necessarily 
arrive  at  the  plane  of  Slit  B  at  a  distance 

Pz  r0 
z  ”  m  c 

from  the  x-y  plane  and  a  distance  r  =  -/x2  +  z2  from  the 


7 

L,  G.  Lang,  "Angular  Correlation  of  Annihilation 
Radiation  from  Solids,"  Carnegie  Institute  of  Technology 
Report,  Sept.  1956. 

8 

R.  E.  Green  and  A.  T.  Stewart,  Phys.  Rev,  98,  486  (1955). 


Figure  10 

Diagram  of  the  Geometry 


Slit  A 


DIAGRAM  OF  THE  GEOMETRY 
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y-axis  where  z  is  a  constant  for  a  particular 
where 


x  = 


Px  r° 
me 


varies  from  x  =  0  to  x 


_  Px  max  r° 
me 


Pz 


and 


(2) 


If  an  annihilation  pair  at  the  origin  have  a  momen¬ 
tum  component  in  the  z  direction  of  P^.  and  if  one  pho¬ 
ton  passes  through  Slit  A  at  the  point  a,b,  the  other 
photon  must  necessarily  arrive  at  the  plane  of  Slit  B 
a  distance  of  z-b  above  the  x-y  plane  and  a  distance  r 
from  the  y-axis  where 

r  =  \/ 'z-  b)  +  ( x  -  a)  (3) 

and  z  is  again  a  constant  for  a  particular  z.  If  Slit  B 
is  a  set  distance  zQ  above  the  x-y  plane,  then  a  coinci¬ 
dence  can  only  be  detected  if  Zo-bo  <  z  -  b  <  Zo  +  bo 
and  if  -  do  <  x  -  a  <  do  .  This  condition  assumes 
constant  efficiency  of  the  scintillators  in  the  x  direc¬ 
tion  and  in  the  z  direction. 

If  the  condition  that  -Oo  <  x  -  a  <  do  is  satis¬ 
fied,  coincidences  will  begin  to  occur  for  z  =  Zo-2bo, 
will  increase  until  z=z0,  and  will  decrease  to  zero 
at  z=zo+2bo  »  Hence,  the  resolution  curve  for  an 
annihilation  at  the  origin  will  have  a  base  width  of  4bo 
centered  about  z  =  zo  .  The  increase  is  identical  to  the 
decrease  because  of  the  symmetry  about  zQ  and  is  linear, 
assuming  constant  efficiency  of  the  scintillators  in  the 
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z  direction.  Therefore,  the  resolution  function  will  be 
an  isosceles  triangle  with  a  base  width  of  4b0  or  twice 
the  width  of  the  slit  in  the  z  direction. 

This  resolution  function  can  be  generalized  for  a 
target  of  finite  thickness  in  the  following  manner. ^ 
Consider  a  target  having  a  width  of  2t  in  the  z  direction. 
A  single  isosceles  triangle  of  unit  height  centered  at 
zQ  is  given  by 


f(Zo,2)  = 


/z-z»/ 

2b0 


G 


jz-Zo|  <  lbD 
iz-Zo\>Zb0  (4) 


so 


}f(2o,Z) 
3  Z  o 


Jf(Zo.Z)  z 


(5) 

(6) 


Hence,  the  total  resolution  function  is  given  by 

R(z),Jf(z.,z)i£  =  J(V 

-t 


2bc 


(7) 


The  function  must  be  integrated  in  three  parts. 
Z+ib0 


lb* 


(8) 


9 

This  follows  the  derivation  given  by  L.  G.  Lang. 
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;  t-Zba<z<-t 


(9) 


-t<Z±t 


(10) 


The  result  of  equations  S',  9  and  10  is  plotted  in 
Figure  11  for  t  =  0.12  milliradians  and  bQ  =  .2  millira- 
dians.  Figure  11  gives  the  total  resolution  curve  in  the 
z  direction  for  a  target  that  is  0.24  milliradians  wide 
and  for  slits  that  are  0.4  milliradians  wide.  The  curve 
is  almost  Gaussian  shaped  and  has  a  half-width  at  half¬ 
maximum  of  about  0.23  milliradians. 

If  the  condition  -a0  <  x- a  <  a o  for  a  point  target 
is  not  fulfilled,  a  correlation  cannot  be  detected  even 
though  the  photon  pair  has  a  pz  that  lies  within  the 
proper  range. 

The  resolution  function  in  the  x  direction  can  be 
calculated  in  exactly  the  same  manner  as  above.  For  a 
point  target  located  at  the  origin,  it  will  have  a  total 
base  width  of  4aQ,  assuming  constant  crystal  efficiency 
in  the  x  direction.  It  is,  therefore,  possible  to  detect 
a  correlated  photon  pair  having  a  p  in  the  proper  range 
if  .  Extending  the  resolution  function  to  a 


Figure  11 
Resolution  Curve 


for  Method  A 
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finite  target  in  the  same  manner  as  above  will  give  a 
base  width  even  greater  than  4aQ. 

If  2a_>>  — — JUgx — ° —  it  can  be  assumed  that  all 

°  me 

photon  pairs  having  a  pz  in  the  proper  range  will  be 
detected  except  when  a  is  close  to  iaQ.  If  2aQ  is  large 
enough,  the  correlations  that  are  not  counted  because  a 

,  p 

is  close  to  -aD  can  be  neglected.  Since  — -- ^qx-  =  10 
milliradians ,  the  condition  for  making  this  approximation 
means  that  the  slit  length  (2a0)  must  be  much  greater 
than  10  milliradians. 

The  first  step  in  the  analysis  of  the  data  is  to 
correct  for  the  z  direction  finite  resolution  of  the 
apparatus.  This  is  done  by  employing  the  parabolic  me¬ 
thod  of  Eckart10  at  each  of  the  data  points.  The  effect 
of  this  correction  is  to  increase  the  height  of  the 
peak  and  to  sharpen  the  corners  of  the  angular  corre¬ 
lation  curves. 

When  Slit  B  is  displaced  an  amount  z  above  the  axis, 
the  corrected  coincidence  counting  rate  D(  pz  )  is  pro¬ 
portional  to  the  total  number  of  pairs  annihilating  per 

unit  time  having  p  =  —  yoZ  .  When  2an  >>  - P*  mox  r°  , 

z  r°  °  me  * 

the  following  relationship  exists: 

oo  in 

D(Pz)  =  const.  J  dppdf  (1 

O  o 


10 


Eckart,  C. ,  Phys.  Rev.  51,  735  (1937). 
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where  ^>(p)  is  the  momentum  space  density  of  annihilating 
pairs  and  the  integral  is  taken  over  the  plane  of  con¬ 
stant  pz  .  It  is  assumed  that  variations  in  p(p)  axe  small 
within  a  region  A  z  corresponding  to  the  slit  width  and 
that  /t>  ( p)  is  spherically  symmetrical  about  the  origin  in 
momentum  space.  Now 


P p  -  \/p2  -  P 


(12) 


so 


<iP  - 

and  a  transformation  of  the  variable  p  to  p  in  the  in- 

r 


(13) 


OO 


tegral  gives 

D(pt)  -  I*  const  \  /o(p)  pdp 

’  A 


since  when  p  =  0  ,  p  =  p  and  when  p  =  OO  ,  p  =  00  . 

r  z  r 


(14) 


Equation  14  can  be  differentiated  with  respect  to  p 


which  gives 

=  -nr  const.  p(pz)Pt 

c4Pi 

and  solving  for  /0(p  ), 

,  ,  1  oLDlh) 

/  =  Z7f  const  pz  (ipL 

Now  pz  =  ^,  so 


(15) 


(16) 


r0x 


L-nm*ck  const 


dD(Pz)  _i 

dz  z 


(17) 
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Since  p(p)  was  assumed  spherically  symmetrical  about  the 
origin,  /0(pz)  =/°(px)  =  /0(py)  '  /O(p)  .  Hence, 


const 


dDiPi)  _l 
dz.  z 


The  momentum  distribution  N  (p)  is  given  by 


(18) 


N  (p)  -  tfrrpl/°Cp) 


(19) 


where  4 IT pl  is  the  surface  area  of  a  sphere  of  radius  p 
in  momentum  space.  Substituting  equation  18  into  equa¬ 
tion  19  gives  the  desired  relation  between  corrected 
data  and  momentum  distribution.  The  preceding  derivation 
is  due  to  Stewart,  ^ 

Z  _  dP(Pi)  z 

const.  dz.  (20) 


When  the  length  of  the  slit  (2aQ)  is  not  much  greater 

than  — Y2~  then  equation  20  does  not  hold  for  the 

me 

relation  between  corrected  data  and  momentum  distribution. 
In  this  case,  the  corrected  counting  rate  D(pz)  is  re¬ 
lated  to  the  total  number  of  pairs  annihilating  per  unit 
time  having  p^  =  by  the  equation 

-  a*, -a  00 


D(Pd  =  (  j 


-a0  -oo 


(21) 


where  a  is  the  distance  from  the  y«z  plane  that  one  of 
the  photons  of  each  detected  annihilation  pair  entered 
Slit  A,  Without  specific  knowledge  of  the  dependence  of 


11 

Stewart,  Ac  T„  ,  op.  cit 
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(°CP)  Px  and  P y  u  is  impossible  to  reduce  equation  21 
any  further. 

B.  Method  B 

This  method  has  been  used  by  several  authors^  *  ^  ^ 
to  analyze  angular  correlation  data  obtained  by  instru¬ 
ments  employing  cylindrical  geometry  but  can  also  be 
used  for  slit  geometry.  The  major  problem  is  obtaining 
the  resolution  of  the  instrument  as  a  function  of  0  and 
Go  where  Go  is  the  correlation  angle  of  the  instrument 
and  7r-^?is  the  angle  between  annihilation  photons. 

Since  Q  and  Go  will  always  be  very  small,  the  fol¬ 
lowing  approximations  can  be  made.  (See  Figure  10 ) 


y  =:  r0  cos  Go  =  ro 

(22) 

Y  -  t'o  s‘n  G-  ==  roG 

(23) 

z0  =  r0  s,n  Go  -  r°G° 

(24) 

If  an  annihilation  pair  at  the  origin  have  correlation 
angle  Q  and  if  one  photon  passes  along  the  y-axis  through 
the  center  of  Slit  A,  the  other  photon  must  necessarily 
arrive  at  the  plane  of  Slit  B  at  a  distance  r  S  Y0  9 


12 

Millett.  W.  E. ,  and  Castillo-Bahena,  R. ,  Phys.  Rev. 
108.  257,  (1957). 

13 

Castillo-Bahena,  R. ,  Ph.D.  dissertation,  The  Uni¬ 
versity  of  Texas,  (1957)  unpublished. 

14 

Deiterman,  L.  H. ,  Jr.,  M.A,  thesis,  The  University 
of  Texas,  (1957)  unpublished. 
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from  the  y-axis  where 

r  =  /x"*+  z1'  -  &-Yb 


n 


(25) 


me  me 

If  Slit  B  is  set  at  a  distance  zQ  above  the  x-y  plane, 
then  a  coincidence  can  be  detected  if  Z0-btf  <r±  /(z4+-b0)V 2/.' 


If  an  annihilation  pair  at  the  origin  have  a  corre¬ 
lation  angle  9  and  if  one  photon  passes  through  the 
opening  of  Slit  A  at  the  point  a,b,  the  other  photon 
must  necessarily  arrive  at  the  plane  of  Slit  B  (y-rQ) 
along  the  perimeter  of  a  circle  of  radius  r-r Q9  and 
with  its  center  located  at  the  point  -a,-b.  If  Slit  B 
is  set  at  a  distance  zQ  above  the  x-y  plane,  then  there 
is  a  chance  for  a  coincidence  being  detected  if  some 
portion  of  the  perimeter  of  the  circle  overlaps  Slit  B. 

Assuming  an  isotropic  momentum  distribution,  the 
probability  that  one  photon  of  an  annihilation  pair  lo¬ 
cated  at  the  origin  and  having  correlation  angle  9  will 
go  through  the  area  Aa4b  at  the  point  a,b  in  Slit  A 
will  be 

A3  Ab 

^TT^oX  (26) 

Let  the  probability  that  if  one  photon  of  an  annihi¬ 
lation  pair  located  at  the  origin  and  having  correlation 
angle  9  goes  through  the  area  AaAb  at  the  point  a,b  in 
Slit  A,  the  other  will  pass  through  Slit  B  when  it  is 
set  at  an  angle  60  be  called  P(B).  Then  the  resolution 
function  R(  6 ,  d0)  of  the  instrument  for  a  point  target 
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located  at  the  origin  when  Slit  B  is  set  at  an  angle  90 
will  be  given  by 


P(B)  will  depend  on  90  ,  9  »  and  the  point  a,b,  at 
which  one  of  the  photons  of  an  annihilation  pair  enters 
Slit  A.  ^  it  wiH  be  assumed  that  for  each  particular 
9  ,  G0,  a,b  that  P(B)  is  equal  to  the  ratio  of  the  arc 
of  the  circle  whose  radius  is  r  -  tq9  that  overlaps 
Slit  B  to  the  total  perimeter  of  the  circle. 

There  are  nine  different  ways  in  which  the  arc  of 
the  circle  can  overlap  the  slit  depending  on  the  values 
of  6  ,  9o  >  a,  b,  and  the  dimensions  of  the  slit.  These 
nine  cases  and  the  calculation  of  P(B)  for  each  are  shown 
in  Figures  12-20. 


15 

This  is  only  true  when  the  annihilation  takes 
place  at  the  origin  where  the  coordinates  are  0,0. 


P<B) 

A 


Figure  12 

»■  No  Arc  Falls  Within  Slit  B 


Figure  13 

P(B)g  -  Arc  Intersects  Slit  B 
Twice  on  Lower  Edge 


Figure  14 

P(B)C  -  Arc  Intersects  Slit  B 
Once  on  Lower  Edge  and 
Once  on  Side 


Figure  15 

P(B)^  -  Arc  Intersects  Slit  B 
Twice  on  Loxver  Edge  and 
Twice  on  Upper  Edge 


p(B)D 
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-r) 
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-V) 
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a 

II 

CO 

♦ 
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P(B)d  . 

<f>r 

2  7 r  r 
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r 
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Figure  16 

P(B)e  -  Arc  Intersects  Slit  B 
Twice  on  Upper  Edge, 

Once  on  Lower  Edge  and 
Once  on  Side 


s 


Figure  17 

P(B)  -  Arc  Intersects  Slit  B 
F 

Once  on  Upper  Edge  and 
Once  on  Lower  Edge 


P(B) 


Figure  18 

-  Arc  Intersects  Slit  B 
Twice  on  Upper  Edge  and 
Once  on  Each  Side 


P(B)S 
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e 
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O 
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II 
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Figure  19 

P(B)  -  Arc  Intersects  Slit  B 

H 

Once  on  Upper  Edge  and 
On  One  Side 


P(B) 


H 


2  7T 


[sin-1  ( 


Qo+a 


)  “  cos’1  ( 


z«  +  b, 


>] 


Figure  20 


No  Arc  Falls  Within  Slit  B 
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For  each  value  of  ,  it  is  necessary  to  divide  9 

into  nine  ranges  in  order  to  integrate  equation  27.  Each 
of  these  ranges  includes  one  or  more  of  the  special  cases 
of  P(B) .  It  is  possible  to  calculate  R(  0  ,  0o  )  using  only 
the  approximations  of  equations  22,  23  and  24  for  six  of 
the  ranges  of  9  .  These  ranges  are  found  to  be  more 
than  sufficient  to  define  the  important  portion  of  the 
resolution  curve. 

In  the  following  equations,  r,  zq,  bQ,  aQ  are  mea¬ 
sured  in  milliradians  so  that  r  corresponds  to  9  and 
zQ  corresponds  to  90  .  Also  let  M  =  So'  \Jrx-  cx>  , 

N  =  3-0 "  yV1-  ,  and  C  =  .  For  \r  -  Z0~Zb0[ 

rb*  ra' 

R(e,e.)  =  c  j  j  p(b)„ 

1  -k.  o  (28) 

For  Z0-  lb.  <  r  <  Z0  : 

«30 

f Uo,90)n  =  c  j  J  clt3cii> 

-b0  0 


+  c. 


P(B)d  dadb 


Zo-K-r 


(29) 


bo 


Z oK-r 


a0 

( 

P(8)c 


M 
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For  Z0  <  r  <  Z0  +  Zbe  ; 


z.+L-r  a< 


R(fl !,«.)„=  cj  jPCBJs^  +  cj  Jp(B)caj<ii> 

—  b.  'o  _  L  /. 

*  °o  M 

b«  ^  fJ 

+  C$  +  C  j  jnB)edaU+c{  (p(B)F  dUb 

Zc^Vr  o  z0+y0-rL  Za+L-r  V 


r  M 

For  z0  +  l\>e  <  r  <  ^Z«»-Ata;*+  <9«a'  .* 


(30) 


“>  ,  A< 


=  c 


V* 


p(8)n  ctddlj  -h  C  P(B)e  cbcU 


> 


L  ? 


+  C 


|  |  P(B)f  dadb 


'L  > 


(31) 


For  +  fxa„;r'  £  r  £  /z/+  Ua.;*-'  ; 

t.  a.  i. 

r69,&)  =c(  jp(B)H&Ji  +cj  pos^d^M 

~K  fe^-a,  -t  7 


o  0 


(32) 


For  f(z0tzbc)i  +■  (I^F  <  r  ; 


b,  <?. 


R(e,6„)E  =  c 


cHi  d  b 


6 


(33) 


63 


and 


There  are  three  ranges  of  r  between  y(z0-Zba)1  +  3c3- 

^(z0'Xbc)x  +  (22»)1  that  contribute  nothing  of 

interest  to  the  R( 9  ,  90 )  curve  and  since  they  are  diffi¬ 
cult  to  evaluate  they  will  be  neglected. 

Since  P(B)^  “  0,  the  first  range  of  r  gives  a  value 
of  zero  for  R(@  ,  90)^.  Similarly,  since  P(8)j.  =  0,  the 
last  range  of  r  gives  a  value  of  zero  for  R( 0  ,  0o  )  .  The 

complete  calculation  of  R(  0 , 0o)jy  is  shown  in  the 
Appendix.  The  same  type  integrals  with  different  limits 
also  had  to  be  evaluated  for  the  other  ranges  of  r.  These 
calculations  give  the  following  expressions  for  R( 6  ,  &o). 
For  Y  ±  Zo-lbt  : 


R(e,e=)i  =  0 


For  Z0 <  r  <  Zt 


(34) 


=  cCC”'e<]  '  (zho~z°  +lr)‘ 


(35) 


For  zrt  £  Y  ±  Z0  +  ; 


~h 


rSc 


xn*r6l 


_  1  _  1  —j 

*  fZYcos  y-z/TTT? 


(36) 
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In  Figure  21,  the  resolution  function  is  plotted 
for  aQ  =  6.35  milliradians ,  bQ  =  .2  milliradians ,  rQ  =  1000 
milliradians ,  and  zQ  =  .6  milliradians  and  3.6  millira¬ 
dians  to  show  how  the  function  varies  for  different  po¬ 
sitions  (z0)  of  the  movable  slit.  These  values  of  the 
parameters  correspond  to  an  apparatus  with  slits  one  inch 
long  by  .016  inches  wide  and  each  slit  a  distance  of  2 
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meters  from  the  target.  The  resolution  function  is  plotted 
in  Figure  22  for  an  apparatus  with  8  inch  long  slits  when 
zc  is  again  .6  milliradians  and  3.6  milliradians.  Figure 
23  shows  the  zQ  =  . 6  milliradian  curves  of  both  Figure  21 
and  Figure  22  plotted  on  the  same  scale  to  compare  areas 
and  Figure  24  shows  the  zQ  =  .6  milliradian  curves  of 
both  Figure  21  and  Figure  22  with  the  maximum  points  nor¬ 
malized  to  compare  shapes. 

No  attempt  was  made  to  generalize  the  resolution 
function  to  include  targets  of  finite  size.  This  can 
probably  be  done  satisfactorily  by  employing  the  methods 
of  numerical  integration  and  will  certainly  broaden  the 
curves.  The  generalized  resolution  curves  for  each  po¬ 
sition  of  the  movable  slit  can  then  be  used  to  correct 
the  raw  data  using  the  general  method  developed  by 
Eckart1^. 

The  corrected  coincidence  counting  rate  at  Q  is 
proportional  to  the  total  number  of  annihilation  pairs 
per  unit  time  having  a  p^>  component  of  momentum.  The 
relation  is 


(°°  (  1  n  l 

D(P/»)  const  )o  J  dfdpy 


where 


Pp=  (W+ft 


(40) 


(41) 
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Eckart,  C,  opf  clt. 


Figure  21 

Resolution  Curves  for  Method  B 
Qq  =  0.6  mrad.  and  3.6  mrad. 
aQ  =  6.35  mrad. 


0.  6  mr. 


(MR.) 


Figure  22 

Resolution  Curves  for  Method  B 
0Q  =  0.6  mrad.  and  3.6  mrad. 
aQ  =  50.5  mrad. 


Figure  23 

Resolution  Curves  for  Method  B 
0O  =  0.6  mrad. 

aQ  =  6.35  mrad.  and  50.5  mrad. 


(Drawn  to  Scale) 
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DC 


CD 
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Figure  24 

Resolution  Curves  for  Method  B 
0O  =  0.6  mrad. 

aQ  =  6.35  mrad.  and  50.5  mrad. 
(Normalized  so  that  the  Maximum  Points  are  Equal) 


74 


Equation  40  can  be  written  as 

oo  an 

PW  =  const  $  (  /°(P )MdPy 

?P  °  >0  (42) 

Since  an  isotropic  momentum  distribution  has  been  assumed, 

equation  42  can  be  integrated  over  to  get 


-  const  zn  \  /O(p)  df>y 
?f>  J0 

(43) 

From  this  relation  Bahena^  has  shown  that  the  distribution 
of  momentum  is  given  by 


Nip) ~  MiPz) 


oo 


o 


dDle) 

oL  9 


where 


ync 


(44) 


(45) 
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CHAPTER  III 


EVALUATION  AND  COMPARISON  OF  METHODS 
A.  Speed  in  Taking  Data 

The  study  of  the  angular  correlation  of  annihilation 
radiation  is  always  a  rather  time-consuming  experiment 
since  the  solid  angles  utilized  by  the  crystal  scintil¬ 
lators  have  to  be  very  small,  of  the  order  of  10'“*’  stera- 
dians.  The  decay  of  the  source,  the  possibility  of  the 
apparatus  getting  knocked  out  of  alignment  while  taking 
data,  and  drift  in  the  electronic  circuits  make  it  ad¬ 
visable  to  collect  the  data  in  as  short  a  time  as  possible. 
This  becomes  even  more  desirable  if  the  target  is  being 
held  at  low  temperature  or  under  a  vacuum. 

The  factor  of  speed  in  taking  data  is  the  big  advan¬ 
tage  of  the  multiple  coincidence  channel  apparatus.  If 
four  properly  spaced  slits  are  used  at  each  end  of  the 
apparatus,  it  is  possible  to  collect  data  for  16  angles 
in  the  same  time  that  a  single  coincidence  channel  appa¬ 
ratus  with  the  same  size  slits  requires  for  one  angle. 

The  data  must  be  taken  with  frequent  change  of  the  angle 
between  the  two  sets  of  slits,  however,  because  the  re¬ 
sponse  of  a  single  channel  can  vary  slightly  due  to 
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electronic  drift  and  temperature  changes,  so  that  the 
sensitivity  cannot  be  considered  as  constant  from  day 
to  day.  By  changing  the  position  of  the  movable  slits 
stepwise  up  and  down  again,  part  of  a  coincidence  curve 
is  obtained  for  each  counter  pair.  Independently  of 
single  channel  sensitivities,  these  16  curves  can  then 
be  joined  together  to  make  one  curve. 

The  usual  arrangement  of  the  photomultiplier  tubes 
and  slits18’19  in  this  type  apparatus  restrict  the  max¬ 
imum  length  of  the  slits  to  about  2  inches.  No  such 
restriction  exists  for  the  single  coincidence  channel 
apparatus;  so  it  is  possible  to  use  slits  much  greater 
than  2  inches  long  and  thereby  reduce  the  counting  time 
advantage  of  the  multiple  coincidence  channel  apparatus. 

B.  Structural  Complexity 

The  SCC  apparatus  (single  coincidence  channel)  is 
easier  to  construct  than  the  MCC  apparatus  (multiple  coin¬ 
cidence  channel).  The  MCC  apparatus  requires  the  con¬ 
struction  of  four  times  as  many  precision  slits  and 
photomultiplier  tube  supports.  The  design  of  the  lead 
shielding  is  necessarily  more  complicated. 

C.  Electronic  Components 

The  MCC  apparatus  requires  far  more  circuitry  than 
the  SCC  apparatus.  It  requires  four  times  as  many  tube 
bases,  preamplifiers,  and  pulse  shapers;  sixteen  times 

•^Trumpy,  Gr  ,  op.  cit. 

19 

Lovseth,  j.,  op.  cit. 
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as  many  coincidence  circuits;  and  a  complex  switch  for 
directing  the  output  of  each  coincidence  circuit  to  the 
correct  sequence  of  registers. 

D.  Expense 

The  MCC  apparatus  costs  more  to  build  than  the  SCC 
apparatus.  It  requires  more  photomultiplier  tubes,  far 
more  electronic  equipment,  and  far  more  man-hours  to  con¬ 
struct. 

E.  Alignment 

The  MCC  apparatus  is  harder  to  align  than  the  SCC 
apparatus.  Eight  slits  have  to  be  adjusted  as  opposed 
to  just  two  for  the  single  coincidence  channel  apparatus. 

F.  Analysis  of  Data 

It  is  much  easier  to  correct  the  data  for  instru¬ 
mental  resolution  by  Method  A  because  the  resolution 
curve  is  symmetrical  and  is  the  same  for  all  positions 
of  the  movable  slit,  while  the  resolution  curve  for 
Method  B  is  asymmetrical  and  varies  greatly  with  the 
position  of  the  movable  slit.  (Compare  Figure  11  with 
Figures  21  and  22.)  It  is  also  much  easier  to  obtain  the 
momentum  distribution  by  Method  A  because  it  only  in¬ 
volves  the  determination  of  slopes  as  opposed  to  the  de¬ 
termination  of  areas  required  by  Method  B. 

In  order  to  Use  Method  A,  however,  it  is  necessary 
that  the  effective  slit  length  be  much  greater  than 
10  milliradians.  As  stated  earlier,  the  maximum  slit 
length  that  can  be  used  in  the  multiple  coincidence 
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channel  apparatus  is  about  2  inches  which  corresponds  to 
about  25  milliradians  at  a  distance  of  2  meters.  In 
order  to  fit  a  standard  photomultiplier  tube  and  still 
cover  the  length  of  the  slit,  the  scintillators  are  usually 
cut  in  the  shape  of  short  right  cylinders.  This  means 
that  the  efficiency  of  the  counter  drops  from  a  maximum 
at  the  center  of  the  slit  to  zero  at  each  edge  which  re¬ 
duces  the  effective  slit  length  by  about  half  and  invali¬ 
dates  the  use  of  equation  21  which  assumes  a  constant 
crystal  efficiency. 

With  a  maximum  effective  slit  length  of  about  12.5 
milliradians,  the  use  of  Method  A  to  analyze  the  data 
obtained  by  a  multiple  coincidence  channel  apparatus  is 
not  entirely  justified.  Several  alternatives  exist: 

1)  use  Method  B  to  evaluate  the  data;  2)  devise  a  way 
to  use  longer  slits  with  the  multiple  coincidence  channel 
apparatus;  3)  collect  the  data  with  a  long  slit  single 
coincidence  channel  apparatus  and  use  either  or  both 
methods  to  analyze  the  data;  and  4)  devise  a  method 
(for  varying  crystal  efficiency)  of  relating  the  count 
rate  to  the  momentum  distribution  from  equation  21. 

At  the  present,  alternative  3  appears  to  be  the 


best  method. 
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The  complete  calculation  of  ft(0>0o)j2  . 
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Each  of  these  integrals  will  be  calculated  separately 
and  the  results  added.  The  integral  for  P(B)p  is 
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Let  these  integrals  be  called  I,  II,  and  III  respectively. 
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Integrating  over  the  variable  a  gives 
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Let  these  integrals  be  numbered  I  through  VII  respectively. 
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Again  there  are  possible  cancellations  in  the  expression 
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